Scaling of Horizontal and Vertical Fixational Eye Movements 
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ABSTRACT 

Eye movements during fixation of a stationary tar- 
get prevent the adaptation of the photoreceptors to 
continuous illumination and inhibit fading of the im- 
age. These random, involuntary, small, movements 
are restricted at long time scales so as to keep the 
target at the center of the field of view. Here we 
use the Detrended Fluctuation Analysis (DFA) in or- 
der to study the properties of fixational eye move- 
ments at different time scales. Results show differ- 
ent scaling behavior between horizontal and vertical 
movements. When the small ballistics movements, i.e. 
micro-saccades, are removed, the scaling exponents in 
both directions become similar. Our findings suggest 
that micro-saccades enhance the persistence at short 
time scales mostly in the horizontal component and 
much less in the vertical component. This difference 
may be due to the need of continuously moving the 
eyes in the horizontal plane, in order to match the 
stereoscopic image for different viewing distance. 



I. INTRODUCTION 

When we view a stationary scene, our eyes perform 
extremely small autonomic movements. These fixational 
(or miniature) eye movements are produced involuntarily 
and are characterized by three different types of move- 
ments: (i) high-frequency small amplitude tremor, (ii) 
slow drift, and (iii) fast microsaccades [1,2]. Generally, 
they serve to counteract retinal adaptation by generat- 
ing small random displacements of the retinal image in 
stationary viewing. Studies of fixational eye movements 
have been going on since the 1950s, but the role of the 
drift, tremor and microsaccadic movements in the visual 
process is not yet fully understood [2,3]. 

When fixating an object its image falls on the fovea, 
the region of highest visual acuity in the center of the 



* jrliang@math.ecnu.edu.cn 
^havlinQophir. ph. biu.ac.il 



visual field. Drifts are slow movements, with a mean am- 
plitude within a range of 1.2 — 9 min arc [2], away from 
a fixation point. Each instance of drift is necessarily ter- 
minated by a microsaccade (cf. Fig.l). Microsaccades 
are rapid small amplitude movements ranging between 
1' and 60' arc and occur at a typical mean rate of 1 to 
2 per second [4]. Microsaccades seem to reposition the 
eye on the target. Tremor (or physiological nystagmus) 
is a high-frequency (ranging from 50 to 100 Hz [2]) os- 
cillations of the eye typically less than 0.01 deg, i.e., less 
than the size of one photoreceptor, and is superimposed 
on drift. Tremor serves to continuously shift the image 
on the retina, thus calling fresh retinal receptors into op- 
eration. If an image is artificially fixed on the retina it 
fades and disappears within a few seconds [5]. Tremor 
causes every point of the retinal image to move approx- 
imately the distance between two adjacent foveal cones 
in 0.1 seconds and thus causes the image of an object 
to constantly stimulate new cells in the fovea [6]. Drift 
and tremor movements are rather irregular and show sta- 
tistical properties of a random walk [7]. Microsaccades, 
however, create more linear movement segments embed- 
ded in the eyes' trajectories during fixational movements. 
There is evidence that microsaccades are (i) persistent 
and anti-persistent at different time scales [3] , (ii) show a 
characteristic signature of suppression and overshoot in 
response to visual change [4,8], and (iii) orient themselves 
according to covert shifts of attention [1]. 

Although finding the specific function of fixational eye 
movements has been a long-standing and controversial 
topic in eye movements research [9], our concern is not 
the purpose of such movements but rather the dynam- 
ical behavior of fixational eye movements and if there 
is some difference between horizontal and vertical fixa- 
tional eye movements. In this article we mainly investi- 
gate these questions using the detrended fluctuation anal- 
ysis [10,11], a technique used to detect possible long-term 
correlations in time series. We find that the persistence of 
horizontal and vertical fixational eye movements exhibit 
pronounced different behavior mostly due to the effect 
of the microsaccades. This result is in good agreement 
with the neurophysiological fact that horizontal and ver- 
tical components of saccades are controlled by different 
brain stem nuclei [12]. Our study indicates that after re- 
moving the microsaccades the scaling behavior of both 
components becomes similar. These findings may further 
elucidate the mechanisms underlying effects of microsac- 
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cades on perception and attention [3,4,8] and their role 
in the neurophysiology of vision [13-16]. In addition, in 
many pathological states the fixation system can be dis- 
rupted by slow drift, nystagmus, or involuntary saccades. 
However, because all three of these occur in healthy in- 
dividuals it may be difficult to determine if there is truly 
an abnormality present. Thus, further characterizing of 
the fixational system may be useful in clinical evaluation 
of such dysfunction. 

II. DATA 

Data was collected from five normal subjects. Eye 
movements for these participants were recorded using an 
EyeLink-II system with a sampling rate of 500 Hz and 
an instrument spatial resolution < 0.005°. The subjects 
were required to fixate a small stimulus with a spatial 
extent of 0.12° or 7.2 arc min (3x3 pixels on a computer 
display, black square on a white background). Each par- 
ticipant performed about 100 trials with a duration of 3 
seconds and total of 474 trials were obtained [3]. The 
recording of each trial includes position trajectories of 
horizontal and vertical components of left eye and right 
eye movements. 
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FIG. 1. Eye position simultaneous recording of horizontal 
and vertical components of left eye movements. The traces 
show microsaccades, drift and tremor in eye position. In the 
horizontal tracing, up represents right and down represents 
left; in the vertical tracing, up represents up and down repre- 
sents down movements. 

Figure 1 shows a typical simultaneous recording of hor- 
izontal and vertical miniature eye movements for the left 
eye from one subject. The horizontal and vertical move- 
ments (upper and lower traces in the figure, respectively) 
exhibit an alternating sequence of slow drift and resetting 
microsaccades. Usually, the subjects show an individual 
preponderance regarding the direction of these drifts and 



resetting microsaccades. In this subject, for example, the 
drift in the horizontal movement occurred typically to the 
right and the microsaccades to the left (Fig.l). 

III. METHODS OF ANALYSIS 

To study the dynamical behavior of fixational eye 
movements we employ the detrended fluctuation analysis 
(DFA) which was developed to quantify long-term power- 
law correlations embedded in a nonstationary time series 
[10]. The DFA method has been successfully applied to 
research fields such as cardiac dynamics [17,11,18-20], hu- 
man gait [21], climate temperature fluctuations [22,23] 
and neural receptors in biological systems [24]. Here we 
apply this method to the velocity series derived from the 
position series of fixational eye movements. For a posi- 
tion series Xi, i = 1, • • • , N + 1, of a horizontal or verti- 
cal movement, we first calculate its velocity series Vi by 
Vi = T(){xi + i — x^, where To is the sampling rate; in our 
experiments To = 500 Hz. We chose to use a two-point 
velocity in order to avoid any smoothing and clearly char- 
acterize the direction and magnitude of a movement. For 
other definitions of velocity see [4]. 

We first calculate the integrated series as a profile 

k 

Y(k) = £>i -<«>], k=l,---,N. (1) 

z=l 

Subtraction of the mean (v) of the whole series is not 
compulsory since it would be eliminated by the detrend- 
ing in the third step [25]. Thus Y(k) in Eq.(l) represents 
actually the "position" . 

We then divide the profile Y(k) of N elements into 
N t = mt(N/t) non-overlapping segments of equal length 
t, where int(iV/i) denotes the maximal integer not larger 
than N/t. Since the length N of the series is often not a 
multiple of the considered time scale t, a short part at the 
end of the profile may remain. In order not to disregard 
this part of the series, the same procedure is repeated 
starting from the opposite end. Therefore, 2N t segments 
are obtained altogether. 

Next, we determine in each segment the best polyno- 
mial fit of the profile and calculate the variance of the 
profile from these best polynomials 

f\ v , t) = - ]T{r(^ - i)t + i) Vv (i)} 2 (2) 

i=l 

for each segment u,v = 1, ■ • • , Nt , and 
1 * 

F>, t) = - Y,{Y{N -(v- N t )t +i) - y v {i)} 2 (3) 

i=i 

for v = N t + 1, • • ■ , 2N t , where y v is the fitting polyno- 
mial in segment v. If this fitting polynomial is linear, 
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then it is the first-order detrended fluctuation analysis 
(DFA1). This eliminates the influence of possible linear 
trends in the profile on scales larger than the segment 
[10]. In general, in the nth order DFA (DFAn), y v is the 
best nth-order polynomial fit of the profile in segment v. 
Therefore, linear, quadratic, cubic, or higher order poly- 
nomials can be used in the fitting procedure. Since the 
detrending of the original time series is done by the sub- 
traction of the polynomial fits from the profile, different 
order DFA differ in their capability of eliminating trends 
of order n — 1 in the series. 

Finally, the fluctuation F(t) over the time windows of 
size t is determined as a root-mean-square of the variance 



crossovers and the exponents at small time scales become 
similar (Figs. 2(c) and (d)). This result indicates that 
microsaccades strongly influence the horizontal compo- 
nents in fixational eye movements. Note, the close sim- 
ilarity of the fluctuation function F(t) in the different 3 
seconds trials, in particular after removing the microsac- 
cades, indicates that the scaling exponent is a stationary 
and significant characteristic of the eye movement. 



F(t) 
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This computation is repeated over all possible interval 
lengths. Of course, in DFA F(t) depends on the DFA 
order n. By construction, F(t) is only defined for t > 
n + 2. For very large scales, for example, for t > A/4, 
F(t) becomes statistically unreliable because the number _ 
of segments N t for the averaging procedure becomes very e 
small. We therefore limit our results to [n, A/4]. 

Typically, F(t) increases with interval length t. We 
determine the scaling behavior of the fluctuations by 
analyzing log-log plots of F(t) versus t. A power law 
F(t) oc t a , where a is a scaling exponent, represents the 
long-range power-law correlation properties of the signal. 
If a — 0.5, the series is uncorrelated (white noise); if 
a < 0.5, the series is anti-correlated; if a > 0.5, the series 
is correlated or persistent. 



IV. ANALYSIS OF FIXATIONAL EYE 
MOVEMENTS 

We applied DFA1-4 to all velocity records derived from 
the horizontal and vertical components. Since the scal- 
ing exponents of the fluctuation functions obtained by 
DFA1-4 are similar, we show here the DFA2 results as a 
representative of the DFA analysis. 

As can be seen from Figs. 2 (a) and (b), the fluctuation 
functions of horizontal components have pronounced dif- 
ferences from the fluctuation function of vertical compo- 
nents. This is expressed by several characteristics, which 
can be observed. There is a broader range of exponents 
in the horizontal compared to the vertical. The crossover 
times, from large exponents (at short time scales) to 
smaller exponents (at large time scales) in horizontal, 
also show a broader range compared with the vertical. 
Moreover, the scaling exponents at short time scales (be- 
tween 12 millisecond and 40 milliseconds) for horizon- 
tal, are typical larger than the corresponding exponents 
for vertical. However, if we remove microsaccades [26] 
the fluctuations of both components, the corresponding 
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FIG. 2. Fluctuation functions obtained by DFA2 for hor- 
izontal and vertical eye movements from the right eye of a 
typical participant: (a) horizontal; (b) vertical; (c) horizontal 
[same data as (a)] after removing microsaccades; (d) vertical 
[same data as (b)] after removing microsaccades. 

In Fig. 3 we show the histograms of the scaling expo- 
nents for the short time scales, for all trials with and 
without microsaccades from the left eyes of all partici- 
pants. From this plot we notice that, at the short time 
scale, the horizontal and vertical components exhibit per- 
sistent behavior (a > 0.5) where the horizontal compo- 
nents are much stronger correlated than the vertical. The 
average value of the scaling exponents for all trials is 
0.76 for the vertical components and 1.1 for the hori- 
zontal components (See Table 1(A)). The scaling expo- 
nents of horizontal components show a broader distribu- 
tion than the vertical components. However, after remov- 
ing microsaccades, the fluctuations of horizontal com- 
ponents and the corresponding scaling exponents have 
a pronounced change to a narrow distribution while the 
vertical components change very little (see Figs. 2 (b) and 
(d); and Figs. 3 (b) and (d)). When comparing the scaling 
exponents for the original horizontal series with the scal- 
ing exponents for the horizontal removed microsaccades 
series, we find that the scaling exponents decrease from 
an average value around 1.1 to 0.74, while for the vertical 
components the scaling exponents decrease from an aver- 
age value around 0.76 to 0.74 (Table I). Horizontal and 
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vertical become similar after removing microsaccades. 
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FIG. 3. Histograms of the scaling exponents a obtained by 
DFA2 at the short time scales [12, 40] ms for all the horizon- 
tal and vertical trials (with duration of 3s) with and without 
microsaccades from the left eyes of all participants: (a) hori- 
zontal; (b) vertical; (c) horizontal without microsaccades; (d) 
vertical without microsaccades. 

We find that at long time scales (between 300 and 
600 milliseconds) , the horizontal and vertical components 
show anti-persistence behaviour (a < 0.5) (see Table I 
(A)), with no significant differences between them. Af- 
ter removing the microsaccades the scaling exponents 
at the long time scales, remain almost the same as be- 
fore. The horizontal components become slightly less 
anti-persistent than vertical (see Table 1(A)). 

TABLE I. Average values of the scaling exponents obtained 
by DFA2 for all fixational eye movements we measured. HL 
= horizontal movements of left eyes, HR = horizontal move- 
ments of right eyes, VL = vertical movements of left eyes, and 
VR = vertical movements of right eyes. 



component 


Short time scale 


Long time scale 




HL 


1.13 ± 0.26 


0.29 ± 0.14 


A 


HR 


1.05 ± 0.25 


0.31 ± 0.14 




VL 


0.76 ± 0.08 


0.34 ± 0.13 




VR 


0.76 ± 0.09 


0.30 ± 0.12 



Microsaccades removed 





HL 


0.74 ± 0.06 


0.26±0.11 


B 


HR 


0.73 ± 0.05 


0.26±0.10 




VL 


0.74 ± 0.05 


0.36±0.14 




VR 


0.74 ± 0.04 


0.35±0.13 



We thus conclude that, microsaccades in the horizon- 
tal components are more dominant than in the vertical 
direction in fixational eye movements. The microsac- 
cades enhance the persistence mostly in the horizontal 
components at the short time scales. At the long time 
scales both horizontal and vertical components are anti- 
persistence and less affected by the microsaccades. 

To further test if the above results are indeed affected 
by microsaccades, we randomly removed parts of the se- 
ries under study with the same length as the removed 
microsaccades and repeated the DFA analysis. We found 
that this procedure does not influence the scaling expo- 
nents. Thus, the scaling difference between the series 
with and without microsaccades is indeed due to mi- 
crosaccades. 

Finally, we tested if the effect of microsaccades can be 
seen also in the power spectral density. To this end we 
analyzed the power spectra of the horizontal and vertical 
velocity series, for the right eye of a typical participant, 
for all trials with and without microsaccades. Results 
are shown in Figs. 4 (a) and (b) where the microsaccades 
are included. The power spectral density of horizontal 
and vertical components are found to be different (Figs. 4 
(a) and (b)). After removing the microsaccades the com- 
ponents become similar (Figs. 4 (c) and (d)). This find- 
ing also indicates that the effect of the microsaccades in 
the horizontal component is stronger than in the verti- 
cal. Note, that this effect is seen much clearer in the 
DFA curves where only a few trials (of 3 sec) are suffi- 
cient to distinguish between the horizontal and vertical 
eye movements. 
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FIG. 4. Power spectral density for the velocity series de- 
rived from the horizontal and vertical components from the 
right eye of one typical participant, (a) horizontal with mi- 
crosaccades; (b) vertical with microsaccades; (c) horizontal 
after removing the microsaccades; (d) vertical after removing 
the microsaccades 
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V. DISCUSSION 

When the visual world is stabilized on the retina, visual 
perception fades as a consequence of neural adaptation. 
During normal vision we continuously move our eyes in- 
voluntarily even as we try to fixate our gaze on a small 
stimulus, preventing retinal stabilization and the associ- 
ated fading of vision [1]. The nature of the neural activity 
correlated with microsaccades at different levels in the vi- 
sual system has been a long standing controversy in eye- 
movements research. Steinman [27] showed that a person 
may select not to make microsaccades, and still be able to 
see the object of interest, whereas Gerrits & Vendrik [28] 
and Clowes [29] found that optimal viewing conditions 
were only obtained when both microsaccades and drifts 
were present. Since microsaccades can be suppressed vol- 
untarily in high acuity observation tasks [30,31], it was 
concluded that microsaccades serve no useful purpose and 
even that they represent an evolutionary puzzle [4,9]. 

Our study using DFA suggests that microsaccades play 
different roles on different time scales in vertical and hor- 
izontal components in the correction of eye movements, 
consistent with [3]. Moreover we show that due to mi- 
crosaccades there is also different scaling behaviour in 
horizontal and vertical fixational eye movements. Our 
results suggest that microsaccades at short time scales, 
enhance the persistence mostly in horizontal movements 
and much less in the vertical movements. 

Our findings that the persistence in horizontal and ver- 
tical fixational eye movements, which are controlled by 
different brain stem nuclei, exhibit pronounced different 
behavior also show that the role of microsaccades in hor- 
izontal movements are more dominant. These findings 
may provide better understanding of the recent neuro- 
physiological findings on the effects of microsaccades on 
visual information processing [13-16]. 

We thank Tomer Kalisky for his assistance in apply- 
ing the DFA method, and Samuel Ron and Lance M. 
Optican for constructive suggestions on the manuscript. 
This work is supported in part by grants KL955/6 
and KL955/9 (Deutsche Forschungsgcmcinschaft) and by 
NSFC (No. 10271031) and the Shanghai Priority Aca- 
demic Discipline. 
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